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SUMMARY 


A numerical analysis was developed to determine the airloads of 
helicopter rotors operating under near -hove ring flight conditions capa- 
ble of producing impulsive noise. The dv-'velopment of the transient 
wake and airload analysis has resulted in a computer program in iidilch 
the solution for the hovering rotor tip vortex geometry, inflow, aero- 
elastic response euid airloads are solved in a coupled manner at sequen- 
tial time steps, with or without the influence of an imposed steady 
ambient wind or transient gust. The computer program was developed for 
future applications i*'* which predicted airloads would be incorporated 
in aji acoustics analysis to attempt to predict euid analyze impulsive 
noise (blade slap). 

The computer analysis was applied to a hovering full-scale rotor 
for which impulsive noise was recorded in the presence of eunblent wind. 
Although blade airloads data were not available for this rotor, the 
results from the initial applications of the conqjuter analysis, as ap_ 
plied to this rotor operating with and without a steady ambient wind, 
are somewhat encouraging. The predicted tip vortex coordinates are in 
reasonable agreement with the limited test data, euid the blade airload 
solutions converged to a periodic behavior for an imposed steady ambient 
wind condition. However, the requirement for refinements to the euialysis 
is evident, particularly for impulsive noise applications where high 
harmonica of blade airloadlng axe required. Recommended refinements 
and suggestions for continued evaluation of the analysis are presented. 
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nn’RODUCTICHi 


The impulsive noise due to blade-vortex interaction, ccounonly known as 
blade slap, is a dominant source of helicopter rotor noise. An understanding 
of its mechanism of generation and a method for predicting it are necessary if a 
reduction in aircraft annoyance and detection are to be achieved. The conner- 
clal and military applications of the helicopter dictate as quiet a machine 
as possible, yet certain operating regimes give rise to impulsive noise levels 
high enovigh to Jeopardize a military mission or be uncomfortable to commercial 
passengers. 

Various investigations have indicated that blade slap arises frcm the 
unsteady lii't fluctuation on a blade due to the induced flow field caused by 
the tip vortex of another blade. An investigation into the noise generation 
of a full-scale hovering rotor subjected to low-speed wind conditions is 
reported in Ref. 1. The results of this test program indice,te that blade slap 
depends on blade-vortex inte^ction emd ccmpresslble conditions at the tip. 

All the test conditions of Ref. 1 had blade-vortex intersections or near inter- 
sections and yet did not necessarily produce iminilsive noise. It was subse- 
quently found (Ref. 2) that the impulsive noise boundary coincided very closely 
to conditions beyond which local blade sonic velocities occurred. It was con- 
cluded that the impulsive noise due to blade vortex interaction of a single 
hovering rotor required a moving vortex to pass close enough to the tip to 
either drive the tip locally supersonic for a short duration giving rise to a 
rapidly forming and disappearing shock wave or to cause a fluctuating motion 
to EUi already existing shock wave. In Ref. 2, it is postulated from the whirl 
tower test results of Ref. 1 that the criteria for rotor impulsive noise are: 
(1) the vortex from a preceding blaule must pass near the blade producing noise, 
and (2) a shock wave must exist (or nearly exist) on the surface of the blade 
where the vortex passes. Thus, the accurate prediction of impulsive noise 
. aquires an ability to accurately predict the following: 

1. the tip vortex position relative to the blade, 

2. the flow velocities at the blade including those induced by the rotor 
wake (particularly the velocities induced by the near tip vortex), 
and 

3. the associated blade airloads including compressibility effects. 

Whereas rotational, noise can be predicted by existing methods, iorpulsive 
noise prediction is difficult at best. The sensitivity of tie induced velocity 
field in hover to the tip vortex position is well known (Ref. 3)* Fluctuations 
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In this velocity field due to variations in the tip vortex position will cause 
existing snock waves to fluctuate or cause local sonic regions to alternately 
form and collapse. Doth situations give rise to a fluctuating pressure field 
which can c.eate unacceptable noise levels. Since the sound pressure level 
is directly dependent on both the instantaneous force and its timewise rate 
of change, it is ai>parent that a requirement exists for an accurate prediction 
of the tune history of rotor airloads in order to produce impulsive noise. 

Vortex induced blade slap has occurred in hover and low speed conditions 
where the tip vortices pass close to the blades producing a major increment to 
the total velocity cmd airloading at the blades. Several analyses have been 
developed by Landgrebe at the Itoited Technologies Research Center (ifTRC) for 
computing the aerodynamic interactions between helicopter rotor blades and 
the rotor wake for the determination of blade airloads and response charac- 
teristics (e.g., Refs. 3 through 5)» A ITTRC method for the hover condition 
is reported in Refs. 3 and 4. This method assumes the hover condition to be 
steady, with the tip vortex pattern from each blade similar and constant with 
time producing steady airloads. It has been hypothesized for hover that a 
small amount of ambient wind is required to produce the unsteady airloading 
required for impulsive noise. The ambient wild may itself be unsteady (a 
transient gust) which results not only in an unsteady variation of the wake 
geometry and airloads but also in a nonperiodic variation requiring a tran- 
sient solution. Techniques used in the above references, were thus applied 
and extended to develop a new aneLlysis applicable to transient conditions. 

Some techniques developed by the author in Ref. 6 were also applied. The 
analysis considers the rotor to be operating initially in a hover mode and 
then, as the gust is applied, the wake is allowed to convect under both the gust 
and its own self induced velocities. The blade inflow velocities and the time 
histories of the blade loadings are then determined. A simplified aeroelastic 
analysis to determine the first flatwise bending mode contribution is included. 


The objectives of this research program were to: 

1. extend the applicability of the existing methods to predict rotor 
wake geometry. Inflow, and airloads for hovering rotors operating 
under conditions producing impulsive noise, 

2. provide the Government with tiic computer iirogram developed so 
that airloads could be predicted and incorporated in a Government 
acoustics analysis to study the imimlsive noise phenomena, suid 

3. demonstrate and perform an initial assessment of the analysis. 


This report contains a description of the computer analysis developed and 
resiilts from initial ai>plications which were performed to demonstrate and 
obtain an initial assessment of the analysis. Results for a hover condition 
without ambient wuid are followed by results for the same condition with 
ambient wind, for \^ich blade slap was encountered as reported in Ref. 1. 
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LIST or SYMBOLS 


Section lift curve elope 

Longitudinal and lateral cyclic pitch control, 
respectively, deg 

Blade chord, m 

Section lift coefficient 

Constants defined in Ref. 7 which are integral functions 
of the blade mass distribution, chordwise center of gravity, 
and flatwise node shapes and their derivatives. 

Flapping' hinge offset , m 

Gravity, m/sec^ 

2 

n-sec^ 

Blade mass per unit span, — ^ — 

Flapwise modal deflection 

Radius vector from vortex filament to point, P, where 
Induced velocity is determined 

Radius of gust front from center of rotor, m 

Blade radius, m 

Differential length of vortex filament 
Normal and inplane blade loads, respectively, n 
Time, sec. 

Resultant velocity at radial station, m/sec 

Velocity at radial station normal to tip path plane, m/sec 

Velocity at radial station parallel to tip path plane, 
m/sec 


Velocity, m/sec 


g 


e 

X, y, z 
a, g 


a 

6 

\ 

r 

H 

h 

6o 

"g 

♦ 

U) 


n 


Induced velocity at point, P, m/sec 
Quit velocity, m/sec 
Flapwlie deflection 

Certeilan coordinate system coo^yonents 

Included angles between ends of straight line vertex 
segment and point at which Induced velocity Is calculated 

Also, blade section angle of attack, deg 

Rotor shaft angle of attack, deg 

Also, blade f > angle 

Flapwlse mode e 

Circulation, m^/see 

Inflow ratio (see Eq. 19) 

Rotor advance ratio 

Blade linear twist angle, deg 

Collective pitch, deg 

Inclination angle of normal to gust front, deg 
Blade azimuth angle, deg 

Natural frequency of 1^^ flatwise mode, rad./®ec*» 
also frequency of transient gust 

Propeller rotational speed, rad. /sec. 
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Subscripts 


Rigid blade axid system at a blade section (Ref. 7 ) 

Final deformed axis system at a blade section, it 'lading 
bending translation and rotation and twist (Ref. 7 ) 

Flapwi&t mode shape indices 

Distance of center of gravity from elastic axis 
Distance of three quarter chord from elastic axis 
Flapwise mode 

Miscellaneous Notations 
Derivative with respect to f 
Derivative with respect to (time) 
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THEORETICAL DEVELOPMENT 


Numerical Model 

The prediction cf airloads and wake geometries of a rotor under the action 
of a tranolunt gust represents one of the most formidable of aerodynamic prob- 
lems. Since such a gust Is usiially of r^hort duratloi and can be time dependent, 
determination of the ensuing loading and vortex wake notion Is an Inherently 
unsteady aerodynamic problem. The fact that the gust lonvects the vortex wake, 
the motion of which affects the loading continuously In time, specifies the 
problem to be one In unsteady aerodynamics. 

In order to model the prcMem of the rotor and its wake interaction with 
the gust, certain assxsnptions are made. The basic assumption Is that the 
rotor aiid Its wake operate Iri an Invlscld medium, a.nd the disturbance they 
create Is Irrotational ; i.e., the problem Is a potential flow problem. The 
blades can be assumed to ae represented by concentrated bound vortices which 
eunounts to neglecting the details of chordwlse loading. This Is the standard 
lifting line assumption. The wake behind the bleuies Is a surface containing a 
distribution of vorticity. This distributed vortlclty can be subdivided into 
two parts: 1) the trailing vortlclty which Is aligned more or less with the 

local velocity, satisfies the Helmholz vortex continuity criterion at the 
blades and provides pressure equalization at the tips; 2) shed vorticity which 
Is laid down Initially paralxel to the trailing edge and Is determined by the 
timewise variations In bound circulation so as to satisfy the Kutta condition 
at the trailing edge. Although shed vorticity can potentledly be of signifi- 
cance in the determination of the timewise variation of blade loading, as will 
be discussed, acciurate vortex modeling of its effect at the blade without 
exceeding practical conqputer limitations remains to be developed. Shed vorti- 
city was thus neglected at this stage of the development of the analysis. 

The loading at the blades is strongly dependent on the induced velocity 
field of the vortex wake which, in turn, depends on its instantaneous posi- 
tion. The instantaneous position of a wake point depends on the time history 
of the velocity field under which it was convected from the insteuit it was 
deposited in the wake until the current time. This velocity field can con- 
sist of a variety of imposed motions (flight speed, gvist, blade rotation, etc.) 
as well as its own self-induced velocity field. While many flight conditions 
permit neglect of wake self-induced velocities as compared to some larger 
imposed motion, in the hover or neex-hover case wake-induced velocities tire a 
primary source of transport and ctuinot be neglected. Thus, a nonlli.ear flow 
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problem is defined which precludea any i;ope of closed form solution, and 
numerical methods are therefore employed. Flgur) 1 is a computerized repre- 
sentation of the wake beneath a hovering rotor .;cn Ref. 3 > The aerody* 
namlc complexity Is evident from the coeplexlty of this wake geometry. 

The numerical treatment considered uses vortex techniques in which the 
continuous vorticlty distribution In the wake Is approximated by several 
trailing vortex filaments divided Into straight seffnents. The blade bound 
vortices are broken Into a discrete number of radial segment.*, over vhlch 
the loading of each Is assumed constant. 

In the physical process a finite lifting surface moving through a 
fluid, trailing vorticlty, w. jh strength equal to the spanwlse bound clrcu- 
latlot. gradient at the Instant of shedding, Is continuously deposited in 
the wake. The numerical counterpart, which analyzes the lifting sxurface at 
discrete time intex*vals, lays down trailing vortices with Initial lengths 
defined by the distance the surface travels during a time interval. Its 
strength is the bound circulation gradient at the beginnins of a time inter- 
val. The possible difference in trailer strengths that occur because of 
time variable loading is taken up by the shed vorticlty (neglected in this 
analysis) deposited at the beginning of each time interval. To magnify 
the conplexlty of the problem, the Individual wake points see different velo- 
city fields (any time or spatially variant imposed motion as well as self- 
induced), and so "lan move different distances luring a given time Intein/'aJL. 
nie vortex segments are free to deform, strotch, euid rotate. 

Under the basic assumptions jvst stated, a model for the rotor wake 
system under the transient gust can be fonnulated. Initially', the rotor and 
its wake are assumed to exist in some "steady state" condition. At some 
Instant, a gust velocity, which can have components in the rotor pleuie and 
normal to it, begins to move across tie rotor wake system (in the generad 
case as a wave front). Points in the wake are then transported by the gust 
velocity, as well as their own self- induced field, as they enter the gust. 
Also, the blade itsjlf may enter and leave the gust, depending on the speed 
of the gust front. Hiese continuously occurring phenomena lead to an 
unsteady aperiodic loading on the rotor blades. 

A numerical api>roximatlun to the above described model does not treat 
the continuously varying problem, but rather yields a solution only at dis- 
crete time intervals. As a result, the following numerical model is used 
to approximately describe the continuous physical process. The rotor and 
its wake are initially assumed tc exist in an initial "steady state" condi- 
tion. That is, an initial circulation distribution, induced velocity 
distribution, and wake geometry are specified. Then, at some instant of time. 
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t, the guat la aasumed to begin moving acroaa the rotor-%rake ayatem. Wake 
points are aasiuned to be convected witit the induced velocity the points 
incur at time, t, as well as the gust velocity If the points lie behind the 
gust front. This transport velocity is assumed constant over the small time 
increment, ^t. Hence, the wake points move the Incremental dlsl'^ce, ^ At, 
where V Is the sum of all imposed velocities Including the gust and the self- 
induced velocities. During the same time interval the blades are allowed to 
advance a rotation distance At ■ OAt, where At 1 b an azimuth angle Increment 
and 0 is the rotational speed. Trailing vortlclty is deposited from the 
blades along the path of travel. At the end of the time Interval, At, a new 
Inflow at the blades and accompanying circulation distribution and loading 
Is detezmilned. The process Is repeated for as many time steps as deemed nec- 
essary to describe the solution to the problem. Figure 2 shows a representa- 
tive comparison between the continuous wake model comprised of a continuous 
sheet of vortlclty and the segmented wake required for computational analysis. 


Wake Convection 

The wake points are convected under the action of their self-induced 
velocity field, the blade-induced velocity and the gust velocity. 

Wake self-induced velocity ana blade-induced velocity . - The self- 
induced velocity field of the vortex wake and the blade -Induced velocity 
field are determined by repeated application of the Blot-Saveurt law, 

_r_ Txdr 

“'p' W u) 

where dv^p is the differential velocity induced at a point in space, P, 
by the vortex segment of strength, T, and differential length, ds, located 
a distance r from the vortex segment. Figure 3 shows this geometry. From 
Ref. 3> foi* a straight line vortex segment of finite length, the induced 
velocity can be expressed in Cartesltu. component form as follows: 

^ ^[(yp-ya) (Zp-^B)"(yp-7B)(?p’7A)] 

''yp* (^p-^a)-(^p-2b) (*p-*b)] 

''*P= ^ [(^p-'^n) (yp-ye) - (^p-^e) (V^a)] 
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« ( 2 ) 
j*(7p-74) (Vp-yt) 

K ■ (f p - 1 /^) (Ip - Ip*! 

AP ‘[(ip-I^)* ♦ (7p - 7a)* (7p - ^a)*]'^ 

BP . [(*p-K,)* +(yp-y,)* *(*P - 1,)*]''* 


Inherent in the use of straight line vortex segments for modeling the 
curved vortex wake are the problems of i) accounting for the curvature of 
the segments and 2) the fact that a singularity in the Biot-Savart law 
exists when the point at which the induced velocity is calculated lies on 
the inducing elemeiit. 

The error due to neglect of vortex curvature increases as the point at 
which the induced velocity is determined gets close to the inducing filament. 
In the vortex wuke the points at which the induced vel )city is calculated 
are endpoints of the Inducing vortex segments. Thus a relatively large 
error can be expected at these points due to neglect of cxurvature. In an 
attempt to accovint for this effect the approach developed for the n 
Geometry Analysis in Ref. 3 is used. Here it is assumed that the » -nents 
adjacent to the point in question can be approximated by circular arcs and 
thoir effects can be obtained by scaling the influence of adjacent straight 
segments according to the equation: 


V. * V, 2KV, 

I Iq ia 


( 3 ) 


where 


« Induced velocity at point, P. 
v^Q a Induced veloci*:y at P, baaed on straight segments. 
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Via •• Velocity Induced by the wake segment that adjoins 
the wake segment adjacent to the wake point. P. 

K • Scaling factor for curvature, a function of vortex 
core radius .tnd azimuth increment (straight line 
segment length). 

The value of K, included as input in the numerical analysis, is given in 
Fig. 4. 

nie sixigularity in the Biot-Savart law is eliminated by assuming the 
vort<)x segment to really have a small but finite cc^e inside of which poten- 
tial flow considerations do not apply. Recognizing t.^'e cooplexity of the 
real flow inside such a region to be beyond the scope of tbs exterior 
potential flow model, it is assumed simply fiab each filament segment in 
the wake conserves circulation per unit length and maintains a soJ.id body 
rotablon within the core region. Hence, within the core 


•?vn« 




2 COS Oj COS /Sf. 

COiC“Tcos7r^ 


( 4 ) 


where Vj^(r) » Velocity induced at r<rj. by Biot-Savart law 

r^ ■ Core radius 

Qf,P,Qfg, Pg » Vortex segment geometry. Fig. 5 

The vortex core radius must be in] ut in the oi'alysis. 

Finally, it must be stated that the wake indiiced velocity is computed 
in an inertial coordinate system shown in Fig. 6. This coordinate system 
identifies the ♦ « 0® (zero azimuth) position of the blades and also serves 
as a reference system for the gust velocity. 

Gust representations. - Various representations for the transient gust 
are considered, nie positive convention for all the models is shown in 
Fig. 7. The least complex gust is a simple step change in the velocity 
impose'^ on all the rotor emd wake points at some initial instant of time 
and held constant in magnitude and direction for all time. 


’ ‘ X'o 

ton e, . 


( 5 ) 
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where 6^ le the angular orientation of the gust In-plane velocity relative 
to the Inertial reference. This particular gust representation has as an 
analogy In flxed-vlng aerodynamics as the classic Wagner gust representation. 

A more complex model consists of the gust moving across the rotor wake 
system as a tra%*ellng wave so that a finite time roust elapse before the 
entire system Is Immersed In the gust. Hence, a point may or may not be 
affected by the gust accordingly as It Is behind or ahead of the gust front. 
Of this model, two possibilities are considered. One has a cofspletely 
xinlform velocity field behind the gust front so that 



0 before g'ist prisaaga 
ivgjj+ jvgy+ kVg^ , after gust passage 


( 6 ) 


and Vgjj, V , v are all constant. This situation Is shown In Fig. 7. To 
determine ^eth§r or not a point lies behind the gust front, let the radius 
vector normal to the front at any Instant of time be given by 


where 


(t) « fg (0) Vg t 

'g, 

7^(0) « 7g at the initial time 


(7) 


Referring to fig. 7» the coL\actlon of points ]ylng on the gust front satis- 
fies the relation 


Sg-fg cos 9q 

yg-rgSiN0g 


= - tan Og 


( 8 ) 


Then, If (xp,yp) < (xg,yg) the point (xp,yp,z ) lies behind the gust front 
and the point Is thus : nflu^oced by the gust velocity 


^=Vg + kVgj Vg (9) 

where Vg “ i'n plane gust velocity, a constant 

=> Normal gust velocity, a constant 

Another degree of complexity to the traveling wave gust can be added 
by assuming a time dei)endent normal gust velocity of sinusoidal form 


Vgj * b sinwt, 


( 10 ) 


where b is a velocity amplitude, uu a given frequency, and tj^ a "retarded 
time" which measures the time since the gust front parsed the point (x ,y , 
2p ) . Hence the gust velocity is represented by ^ 


_ r 0 (>lp»yp)>(*q'y9) 

'^9 ^ Vg k bsin (ut| (Xp»yp) <(*qiyq) 


(U) 


where are the in-plane Cartesian coordinates of the point in question 

(^ »yg ) in-plane coordinates of points on the gust fr.mt. The 

time, tj^, that the point has spent in the gust, relative to the gust front 
location is 

. l^gl-lfplcOS (fip-®g) 

•i = ' ' ■; — : /in\ 


Thus, a point, is determined to be behind or ahead of the gust 

front, and its velocity is then determined by the position relative to tlie 
gust front, as determined by tj^. 

If desired, the analysis could easily be expanded to include other 
types of gusts such as a ramp gust in which the velocity normal to the tip 
path pleine increases linearly with position in the gust. Hence, the gust 
velocity would be given by 


(0 (>ip.yp)>l>‘q.yq) 

1 Vg + k m 

^ cos(0p-eg)‘ 

[ Vg 


(Xp«yp)<(*9»yq) 


( 13 ) 


where m is the gust acceleration in the direction normal to the tip pach 
plane . 


Thus, in accordance with the gust model desired, wake points convect under 
the gust velocity as well as their own self-induced velocities. 
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Determination of Circulation and Inflov 


In order to determine the wake induced velocity field at 'ny inotemt 
of time, the wake circulation distribution must be known. To determine the 
circulation within the framework of lifting line theory, the fundamental 
clrc'ilatlon eqxiatlon (Eq. l4) is used. This equation is an implicit state- 
ment of the Welsslnger lifting surface theory boundary condition. 


and 


rM)' 


u(r,'l') 


c(r) Cjj(r,'#') 



c(r)a(,(r,'^)a(r,'^) 


(14) 


u(r,«- [u*(r,<-) tulfr.*)]''* 


(15) 


where U(r,t) is the resultant velocity at the radial station, r, of a blade 
located at azimuth location ; c la the local chord; a^ is the lift curve 
slope of the section; a is the local angle of attack. r(r,i^) is the local 
bound circulation at the "r" radleil position on the blade at azimuth angle i/ . 
The circulation is stall limited with a maxlmxim value determined from the 
section steady state stall. The angle of attack (Fig. 8) contains the blade 
pitch angle, any cyclic angles, aeroelastlc motion and Induced velocities, 
due to both the wake and other blades. Ac a result 

a(r.V') »e + tori' Up/Uy 


where 




♦ _ , 

COS V'- Wg COS0 + 

cose-w, sir d\ - rP* 


+ Wg cos 9 



S'ficos f Wj sin 5 


(17) 
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^ Sind 

^r 

- w, d cos d - ^ cos d sin dj + w,' sin d cos ♦ 
” " ^0,- CO* ^s cos d cos*^ 


(18) 


and 




vsin a,-v,j 

“ flR 


(19) 


Equations 17 and 18 are developed fully in the aeroelastic analysis 
of Ref. 7 where they Include flatwise, edgewise auid torsion effects. Here 
they have been modified to Include only flapping emd flatwise bending, and 
lead-lag effects have been neglected. Equation 17 was modified to Include 
the vertical gust component 

. equation it + ^ (20, 


Even though the In -plane gust component should be snail relative to the 
rotational velocities, and so coxild be neglected without serlc.ua effects 
In Uj/flR, It was also Included. Thus 



» Equation I8 + 


Vg»Sin '^'-VgyCOSV' 

HR 


( 21 ) 


The angle, 0, includes blade collective pitch, twist and control system 
Input so that 


0 * do 

dg* d,(5+r -0.75) 

‘«t 5 r"^s cos'^-B.s sin ’I' 


( 22 ) 

(23) 

(24) 
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where 6.75 is the collective pitch at 0.75 radius, linear twist 

over the radius, and and are control itiputs. Pitch-flap coupling 
effects ore neglected. Ihe quantity, Xg, is the inflow ratio and contains 
the induced aerodynamics. As a result it contains induced velocity contri- 
butions fro-.n all bound and wake vortex elements. Equation l4 con be written 
for a number of radial stations on each blade. The aerodynamics of the blade 
and wake system enters through the section angle of attack, a, which contains 
the induced velocity of the WL.ice as well as tlie effect of other blades. The 
system of equations developed by applying Eq. l4 in this manner can be ex- 
pressed as a set of simultaneous equations wltn the unknown boxuid circulations, 
Tj's, at the various radial stations to be determined. 



(25) 


Here the subscript i identifies t\e particular load point at which the 
velocity is induced, J identifies an inducing element of unknown strength, 
Pj. Aj^j is the influence coefficient which incorporates the effect of rela- 
tive geometi7 between the unknown vortex segment and the 1^4 load point. 
Bi represents the known inflow at the i'^4 load point and includes all known 
velocities normal to the tip path plane. This is a fairly standard approach 
and has been used with success for a number of years (e.g.. Refs. 3-5). 


Blade Loads and Aeroelastic Effects 

The aeroelastic effects are developed in detail in Ref. 7* The per- 
tinent equations and the implementation in the computer program are 
discussed here. A major assumption in accovintlng for the aeroelastic 
effects is that flatwise motions and deflections are the major contribution 
for the hover conditions so that edgewise and torsion characterist j cs can 
be neglected. Hence the equations of flatwise bending for a given blade as 
given in Ref. 7 are reproduced here, neglecting the edgewise and torsion 
effects. 
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0 •/ Vwj ®] dF - C,j{q* 4 . Qwl*^*.* 

~Bq~Z$ «o ■*■ /3 •'" 2 flo - "i- (I ♦ C08 2 ^)U + T q^ , 

i ' 

• [-2 e, uJ^* C,^ ., 4. 0, sin 2 flo Ce,^ - 2 ^ sin (C«a* Ce9. j] 
-C, 0 | (« ^ cos 00 •»■ cos 00 - 5y^ sin Q^) - C, 2 ^ [(S +/3) cos 
^Z$0 sin 0o] + C, 5 . (^ *+ /3) sin 00 - Cea^ [i*+ /? 0 “ cos 2 0o) 


4 - sin 2 0 o] - 0 , cos 2 BqIc^q -OH 0 * 3 ^) - 0 , C*Oj 

vrtiere i, i’ are Indices on mode shapes. The C's are constants defined xn 
Ref. 7 which are integral functions of the blade mass distribution, chord- 
wise center of gravity, and flatwise mode shapes and their derivatives. 
Note (hat blade lead-lag motion, pushrod effects and damping have been 
neglected. 


Likewise, for each blade the flap angle equation of Ref. 7 can be 
written, usina the aaine esaianntionfi as the flatwise bending equation. 



r d f - Rg 7og (§ /3 + g - ig / 3 ) - C 97 (^*cos 6 ^ 


where 


- 0 * 


Sin 0„ 


* 

' 2i90o sin 6 ^j- sin 0 o) 



( 27 ) 


Equations 26 and 27 define a set of simultemeous equations from which the 
modal accelerations, and flapping acceleration, V, can be obtained. 

Straightforward nxanerical Integration will yield the modal and flapping 
amplitudes so that the aeroelastic flatwise deflections are given by 


'^e = ? Qw (t.n) y^(r) 
1 1 

( 28 ) 


( 29 ) 


(30) 
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_ A 

wticre w'q la the slope of the deflection and Wq is the velocity due to aero- 
elastic motion, n is an index defiring the particular blade in question, 

I, 2, «•« , B. 

The ijiipiementation of the aeroelastic considerationa represents an 
attempt at close coupling the aeroelastic procedures wltti the aerodynamic 
effects. Since the problem is solved sequentially in time, tlte aeroelastic 
effects and blade loads must also be c'eterriined at each time step, jt should 
be noted that Eqs. 26 and 27 for the modal accelerations and blade flapping 
acceleration at any instant depend on the current value of the modal ampli- 
tudes and flap angles vhich are not known a priori. A sii|iple technique is 
considered here. It is assumed the modal quantities, 6, R, etc., are known 
at time t » C. Based on initial circulations and induced velocities, initial 
modal accelerations are determined. A simple Euler integration scheme (actu- 
ally the kinematic relations for velocity and displacement) serve to yield 
these modal qurntities for the next time step. Then, following blade advance 
and wake def'^rr atlon, new blade loads are determined giving new modal accel- 
erations and the process is repeated for the next time step. 

The blade loads, as developed in Ref. 7 are given by 


“ T ^ i^f cos <#> ■^C^, sin <P) (31) 

i (-^)f (^) (C| sin^ -C^COS<^) (32) 

These can be converted to conventional distributed loads by multiplying 
through by ny)^. 
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{Xirther Assumptions 


Asiuraptions have been made in order that the computer program be feasible, 
both economically and practically in terms of required computer storage and 
run times. Since that part of the analysis which becomes more expensive in 
both the above categories is the itep-by-step analysis of the distorting 
wake, most of the assumptions made sure in this area. 

A primary assumption regarding the distorting wake is neglect of the 
shed vorticity. This amounts to directly neglecting instantaneous increments 
in the bound circulation and thus the aerodynamlu load. The missing incre- 
ments in bound circulation will also be reflected as missing increments in 
wake element strength, the effects of which are treuismitted to the blades 
through the wake induced velocity field. Some Justification for neglecting 
shed vortex elements arises from the fact that the circulation strength of 
the shed vortex elements are normally much lower than the strength of the 
trailing vortex elements. However, in determining an instemteuaeous load 
accurately, as is necessary for determining a high number of loading har- 
monics, neglecting this shed vorticity is potentially serious. Also, the 
presence of shed vorticity can alter the shape of the chordwise loading. 

While this might not be too significant in tenns of average loading, it 
should be important in determining the higher harmonics which are important 
in teimis of the noise signature. Inclusion of shed vorticity should be con- 
sidered for future refinement of the analysis. 

To further sinqplify the problem, it is assumed that the induced effects 
are most described by the tip vortices from the individual blades. Hence, 
the wake distortion is primarily described by the interaction of the various 
tip vortices. The contribution of the inboeurd vortex sheet to the tip vortex 
distorting velocities is assumed constant which implies that the relative 
positioning of sheet and tip vortex does not change much. This permits the 
contribution of the inboard sheets to the tip vortex distortions to be com- 
puted once, initially, euid held constant for all times. Another, perhaps 
more questionable, assumption is that the tip vortex has no effect on the 
Inboard sheet, although a provision for the sheet to convect under the gxxst 
velocity does exist. Essentially what this does if remove the problem from 
the realm of the full interaction problem, for if the sheet were to deform 
under the tip vortex, it would also have to be allowed to deform under its 
own self-induced field vriiich, computationally, is what is being avoided to 
prevent prohibitive computer time requirements . 

The assumption of the wal« vortex segments having a core of finite 
radius with the internal fluid performing a solid body rotation has been 
disc\issed in conjvinction with the development of Eq. 4. It can be restated 
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here, however, that such an assumption is necessary to eliminate a singular* 
ity in the Biot-Savart lav even though the proper determination of this 
quantity is open to conjecture. The vortex core radius is an input quantity 
in the analysis. 

The aeroelastic phenomena are assumed to be adequately modeled by the 
rigid body flapping and the first flatwise elastic modes. For the peurticular 
case of a transient gust passing over a hovering rotor and wake system, higher 
modes and cc tpled effects are most likely negligible. Any errors to higher 
harmonics of r le loading should be less than those provided by neglect of 
shed vorticlt,. and the lifting line assumption. 

Another feature considered to minimize the computational requirements 
is variid)le spanwise spacing on ttie rotor blades depending on the particular 
item of the solution being considered. The radial segments for determining 
wake distortions are not necessarily the sane as those for determining circu- 
lation and inflow. No definite relationship between the number of panels 
for distortion and the number for circulations have been established; the 
final distribution selected was based on obtaining a good representation of 
the tip region. One restriction is that the set of filaments used to compute 
the distortions is a subset of those used to compute the circulations. The 
required circulations for the distortions are then sirrq)le averages over the 
number of ;>anels between successive filaments used for the distortion. This 
amovints to maintaining the same average circulation (or load) between the 
same filaments in both the distortion solution and the circulation solution. 


Sequence of Operations in Computer Program 

The preceding analysis was programmed for application on a UNIVAC 1110 
digital computer. This consisted of exi>anding two existing techniques used 
in UTRC computer progreuns to perform the necessary computations. The two 
parent programs ai^ the Circulation Program and the Distorted Wake Ceoroetry 
Program described by Landgrebe in Refs. 3 » ^ & 5 . 

The Circulation Program computes the circulation distribution over the 
rotor, given the blade control and response characteristics and a specified 
wadee geometry. In its simplest form, the wake is a classic skewed helicoid 
in which the trailing vortex filaments are convected relative to the blade 
by the resultant velocity equal to the vector sura of blade rotational velocity, 
rotor forward velocity and momentum induced velocity. If wake distortions 
are important, the basic Circulation Program can be coipled in ar Iterative 
procedure with the Distorted Wake Geometiry Program. The Distort j. Wake 
Geometry Program is an analytical method whereby, given a bound circulation 


21 


distribution, the distortion of the wake filaments can be confuted. The 
basic approach of this method is straightforward and can be sketched as 
follows: First, an initial wake geometry is assumed and the distribution 

of circulation along each filament is known from the assumed bound circula- 
tion distribution. The classical Biot-Savart law is applied to numerous 
wake points to determine the Induced velocities at these points. These dis- 
torting velocities eure then numerically integrated over a small time increment 
by the iiuler technique to determine new wake element positions. The process 
is repeated until the wake points converge to some defined configuration for 
the given circulation distribution. This wake geometry can then be supplied 
to the Circulation I^ogram and new circulations computed. If the new values 
do not aigree with the old values within a prescribed tolerance, the new cir- 
culations can be supplied to the Distorted Wake Gieometry Program to determine 
a new wake geometry. This wake can then be supplied to the Circulation Pro- 
gram to conqpute still another set of circulations. This process can be 
repeated until two successive sets of circulations converge to within a 
prescribed tolerance. 

The previously described method, by virtue of the way in which the wake 
and circulations are treated, has within it implicit assumptions of steady' 
state or periodic coi.dltions. In general, the treuisient gust problem is an 
unsteady, aperiodic phenomenon so that direct application of this procedure 
cannot model the problem. Thus, for the analysis described herein, the 
physically more realistic but numerically more difficult procedure was used 
which consisted of satisfying flow conditions at each point on the blade- 
wake mesh at each Instant of time. Hence, the distorted wake geometry is 
computed at each time step and the circulations are determined at the indi- 
vidual blade load points at each time step as well. The aeroelastic contri- 
butions to the circulation and inflow are also computed at blade load points 
at each time step. With this fully coupled solution and the creation of new 
wake points at each time step, the core storage and computation time require- 
ments are increased substantially. Following the determination of the 
instantaneous blade inflow velocities and the aeroelastic contribution to 
blade position and velocities, the radial distribution of blade loading is 
determined at each azimuth position (time step). After the rotor has been 
allowi*d to move through its prescribed number of time steps, the blade load 
distribution is Fovurier analyzed in time (azimuth). 
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This procedure is summarized in the foUowin({ steps: 


Procedure >‘or Tims 0 -- Initial Blade Positions 

1. Calculate Blade Coordinates. 

2. Calcxilate Wake Coordinates Vor Starting Hover Condition 
based on Input Wake Constants (Prescribed Cieneralized 
Wake of Ref. 3)« 

3> Calciilate Or Input Circulations And Induced Velocities. 
4. Calculate Blade Response And Airloads. 


I>rocedure For 1-^ch Time Step 

1. Advance Blad&s One Azimuth Increment. 

2. Add New Trailing Wake Elements. 

3. Compute Distorting Velocities And Oust Velocity At Each 
Tip Vortex Point. 

4. Compute Oust Velocity At Each Inboard Wake Point. 

5. Compute Wake Geometry. 

6. Solve For Circulations At Blades Including Effect Of 
Gust Velocity And Blade Response. 

7. Compute Induced Velocitie-J At Blades. 

8 . Compute Total Blade Velocity Components Including 
Gvist Velocities. 

9. Compute Blade Response And Airloads 
(Flapping & 1st Bending Mode). 

In svimnary, a numerical method capable of performing an airloads 
analysis on a hovering rotor subjected to a transient gxist has been developed 
The method computes wake distortion, blade inflow, aeroelastic response and 
airloads within the restrictions in 5 >osed by the assumption described herein. 


DISCUSSION OF RESULTS 


To perform a demonstration and Initial evalxiatlon of the analysis, it 
was applied to a full-scale rotor test condition. The test rotor is a stan- 
dard 3-bladed Boelng-Vertol CH-U7B rotor with a blade radius of 9,lUU m 
(30 ft), a blade chord of 0.61*1 m (25.25 In.) and a blade llnaar twist of 

-9.12 dog. The test was conducted at the Boelng-Vertol Whirl Tower facility. 
The selected test condition is defined completely in Ref. 1 (Run 116). The 
rotor was operated at a nominal thrust of 1.20 mn (27,000 lb) and a tip speed 
of 228 m/sec (7U7 fps). The test condition (Run ll6) produced blade slap in 
the presence of an ambient wind speed of 4m/sec (9 mph). 

The analysis was mn to simulate both am ideal steady hover condition 
(no eunbient wind) autd the ambient wind condition. The ideal hover condition 
was run primarily to (1) check the analysis for convergence by demonstra- 
ting that the solutions (wake geometry, blade response im'i airloauls) for 
hover eure time independent, atnd (2) to evaluate the ac;ur''.y of the waike 
geometry portion of the analysis by comparing the predicted tip vortex geom- 
etry with available experimental information obtained from generalized model 
rotor hover test data. The ambient wind condition was run to (l) predict 
the tip vortex geometry and airloads for an experimental near hover condition 
which produced blade slap, (2) check the convergence of the solutions to 
periodic variations expected for a steady wind condition, and (3) compare the 
periodic time history of the predicted blade-vortex separation distance with 
test data. Unfortunately, measured blade airloeid data was not available for 
comparison. 


Hover Condition 

Initial computer runs with the transient airloaul analysis were made for 
the hover mode in order to test the fully coupled time step method of solu- 
tion. The conditions of Run ll6 of Ref. 1, without the ambient wind, and a 
time step increment defined by an azimuth increment of 30 deg were used in 
the analysis. To start the problem, an initial wake geometry was prescribed 
using the generalized hovering wake coordinates from Ref. 3 for the appro- 
priate number of blades, twist and thrust level. The intent was to deter- 
mine whether or not the calculated tip vortex geometry would return to the 
Initially prescribed geometry following a hopefully converged iterative solu- 
tion of the associated blade and wake circulations, blade response and air- 
loads . 

Since the rotor was assumed to be in a steady hover condition, temporal 
deviations in the botind circulations should not exist in the final solution. 
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However, in the initial application of the analysis, it was found that 
unreasonable oscillations in the value of circulation at the three outboard 
radial stations occurred by the third tine step. No tendency was found for 
these oscillations to diminish as the free wake developed. In investigating 
the cause of these osclllaticns, it was found that the air loading in the tip 
region was very sensitive to the representation of the wsike in the tip region. 
Special attention to the modeling of the wake at the blade tip was thus 
necessary. 

The wake representation, as originally conceived, wsis assumed to have 
the trailing vortex segments positioned ais straight line segments attached 
directly to the bound vortex (Fig. 9)« An unreasonable magnitude of the 
contribution of the trailing vertex element attached to the blade tip prompted 
a reevaluation of this blade-wake representation. The various blade-wake 
representations considered are shown in Fig. 9 in which P (t) implies a cir- 
culation value at the end of a given time step and P (t - implies a value 
at the etid of the previous time step. The various vortex models considered 
represent attempts to improve the modeling of the tip region without going 
to a complete lifting surface model with shed vorticity. The first permuta- 
tion on the original model, representation (2) in Fig. 9» vkb a simple approx- 
imation to estimate the effects of chordwise blade geometry. This amounted 
to assuming the bound vortex to be located at the locail one-quarter chord 
line and modeling the chordwise and trailing vorticity as follows. Vorti- 
city was fixed to the blade in the chordwise direction up to a distance one- 
quarter chord aft of the trailing edge. At this distance, trailing vortex 
elements were initiated with strengths determined by the bound circulation 
gradients at the beginning of the time step. This representation was defined 
by considerations of numerical lifting surface theory. Although this repre- 
sentation was expected to improve the results, the erratic oscillations of 
the circulation and induced velocity solutions persisted. Blade-wake repre- 
sentation (3) in Fig. 9 wB’B then used. This representation is similar to 
representation (2) except that the chordwise vortex elements were omitted. 

On the basis of this investigation, representation (3) was selected as 
the tip model which most reasonably maintained the temporal constancy of 
bound circulation and produced a reasonably accurate tip vortex geometry for 
the hover condition. The convergence of the calculated hover waJie is shown 
in Fig. 10. The experimentally determined hovering wake coordinates from 
Ref. 3 compared with the predicted coordinates after 4l, 62, and 8l time 
steps corresponding to 1.7, 2.6 and 3*4 blade revolutions, respectively. The 
convergen-^e of the near wake < 480 deg) and the good Eigreement of the 
results from the 8lst time step with the experimental model data is evident. 
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particularly in the important region around 120 deg where the following blade 
is located. The far wake U80 deg) is not converged to a steady geom- 

etry. However, this is consistent with the experimental findings rei>orted in 
Ref. 3 which indicate the tip vortex for a hovering 3*bladed rotor is 
stable geometrically beyond approximately UdO deg. 

It should be mentioned that representation (3) does give an upwash field 
at the tip which tends to drive the tip into stall. This is due primarily 
to the reduction of the effect of that part of the self -generated tip fila- 
ment closest to the bound vortex so that the contributi<x\ to the induced 
velocity field from this filament is more on the order of the tip filament 
from the first preceding blade. The resulting upwash is then due primarily 
to the contraction of the tip filament from the preceding blade as it passes 
near the radial load point. 


Hover With Ambient Wind Condition 

As mentioned previously, the rotor operating condition from Run ll6 of 
Ref. 1, including the presence of the wind, was selected as the hover with 
ambient wind condition. The analysis was started using the same prescribed 
experimental wake geometry as the previous ca. e for the hc/er without wind 
condition. The 4 m/sec (9 mph) ambient wind measured for Run 116 was applied 
over the entire rotor at the first time step as an Impulsive gust from the 
225 deg ("northwest" in Ref. 1) azimuth direction. In order to simulate 
suialyt ically the imposed steady wind condition, the rotor was subjected to 
the gust representation described by Eq. 5 with Vgj, ■ 0. The results of 
applying the analysis to this condition are £ nmarized in Figs. 11 through l4. 

Initially calculated time histories of blade airloading (thrust loading*) 
at selected radial stations, 0.80 R, 0.87 R, 0.939 R O.985 R, are pre- 
sented in Fig. 11. In this figure, the blade azimuth angle is representative 
of time; i.e., each 15 deg of blade travel represents a time step in the com- 
puter analysis. The analysis was run for 74 time steps which resulted iri 1110 
deg of azimuth travel for each blade which is slightly more than three revo- 
lutions. As shown in Fig. 11, the predicted blade transient loading tends 
to diminish after about two revolutions of blade travel at all radial sta- 
tions. The maximum values attained at 0.939 R as well as the high, nearly 
constant, values at O.985 R are due to 


♦Blade thrust loading is presented in metric units, N/cm. To convert to 
English units, Ib/ir, multiply by 11. 3. 
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these stations reaching their stall linlt under the influence of the upwash 
induced priinarily by the tip vortex of the preceding bln.de. Although experi- 
mental airloads were not available for compariscn, these initial results 
tended to be as anticipated in that they show an Initial transient response 
followed by a tenlency to settle to a periodic, one-per-rev variation. A 
periodic distribution, similar for each blade, is realistic for thi'' condi- 
tion where the ambient wind, initially applied as an impulsive gust velocity, 
is continued in time as a steady wind. However, the solution was luestioned, 
due to the high harmonics ("spikes”) in the thrust loading (e.g., at 0.872 R 
near ^ ■ 700 deg for blade 3 and at 0.939 R near 780 deg and 1130 deg for 
blade 2). These spikes were found to be caused b/ numerical rather than 
physical meeuis, and were due to an inconsistency in the vortex model, ’^hls 
ln»’onsistency occurred because of the attempt to model the inflow at the 9 
radial stations and interpolate to 1^ radial stations to compute the blade 
loads. Since the 9 radial positions at which the inflow is computed are not 
a subset of the 15 positions at which the blade loads are computed, the inter- 
polation scheme to compute the inflow at the loaui points, given the inflow 
points, can lead to completely ei roneous results. These spikes were traced 
to this inconsistency and so must be discounted as a possible noise source. 

The case was recomputed with the aerodynamic load stations selected to 
coincide with the ') induced velocity stations to provide a consistent model. 
Airload data at four outboard blade stations (0.75 R» 0.85 R» 0.985 R and 
0,99 R) '’or 3 revolutions are presented in Fig. 12. Similar to the previous 
results, the blade loads very quickly settle into a one-per-rev periodic 
fluctuation compatible with a steady forward flight response. The in-plane 
loads, presented in Fig. 13 » show very little unsteady response. 

Comparing Figs. 11 and 12, most of the numerical transients have been 
eliminated by requiring the aerodynamic load points and inflow stations to be 
coincident. It can be noted that the small spike near <|[ « 500 deg in Fig. 12 
is due to a near vortex passage, and may be considered a transient nonrecur- 
ring effect for this steady wind condition. These results show that in 
modeling the lime histories of rotor aii loads by a transient analysis with a 
deforming wake that it is necessarv fnr airload and Inflow stations to b« 
coincident. Even though the comparison does demonstrate the sensitivity of 
the results to the model selected, it must be emphasized that inferences from 
combining the two sets of data must be made with care. The sets of airload 
stations are different, and the 9-station set is not a subset of the 15 orig- 
inal stations. 
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It is recognl;:ed that the selected steady rotor operating condition does 
not represent a true transient gust case. A transient gust Implies a time- 
wise variation in wind over a rotor-wake system In equilibrium. What has 
actually been Investigated Is the condition of a uniform in-plane wind applied 
Impulsively over a rotor-wake system which Is not initially In equilibrium. 

As a result, the transients noted are not realistic, but are numerical due to 
this nonequilibrium. In future applications, when a true gust condition Is 
to be Investigated, the procedure should be to run the analysis in the steady 
condition (without the gust) for a sufficient number of time steps to achieve 
equilibrium of the solutuon, euid then the gust velocities should be applied. 
This should eliminate numerical transients in the predicted transient air- 
loads produced by the gust. 

Due to the limitations in the available experimental data for correla- 
tion purposes, the achievement of the one-per-rev periodicity and the real- 
istic wake results predicted for the hover (no wind) condition have been the 
major items presented to this point which indicate some physical validity of 
the new analysis. Other than acoustic data, the one parameter for which use- 
ful test data is available in Ref. 1 is the time history of the blade-vortex 
separation distance. This is the vertical distance between a blade and the 
tip vortex from the preceding blade as shown in Fig. ''U, and was measured 
using a rotating camera mounted at the hub of the CH-47B rotor. Considering 
the limited extent of the whirlstand data and the measured accuracy of the 
wind conditions, good agreement is shown in Fig. l4 between the analytical 
results ar I the test data. The analytical separation results exhibit a one- 
per-rev sinusoidal variation with blade azimuth similar to that indicated by 
a composite of test data form several runs in Ref. 1. The repeatability of 
the analytical results for blades 1, 2, and 3 is consistent with the repeat- 
ability of the test data foi two (2 ) rotor cycles which is indicative of a 
steady wind condition. The differences in amplitude and phase angle relative 
to the upwind direction are presumed to be at least partially attributable to 
the inaccuracies in the measurements of wind velocity and direction. The 
wind measurements were made on a tower 2.5 diameters (4.57 m, 150 ft) away 
'"rom the rotor center, and the wind direction was recorded only generally as 
northwest for this condition rather than as a specific azimuth in the fixed 
reference system (such as the 225 deg used in the analysis). These results, 
coupled with the correlation demonstrated in Fig. 10 for the hover condition, 
provide eui encoureiging inuication of the capability of the analysis to pre- 
dict reasonably accurate tip vortex geometries. 
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CONCLUDING REMARKD ANT tlBCWf4ENDATIONG 


The development of the transient wake and airload ajialysls has resulte<l 
In a computer program in which the solution for the hovering rotor wake, 
inflow, aeroelastic response and airloads are solved in a coupled manner at 
sequential time steps, with or without the influence of an imposed gust. 

Based on the Initial applications of the analysis, reported herein, the 
predicted tip vortex coordinates are in reasonable f.greement with the limitei 
experimental data. Unfortunately airload data was not available for comperi- 
son, so that only the achievement of the periodic behavior of the airload 
solutions for a steady ambient wind condition could be used to indicate some 
degree of physical validity of the analysis. 

Although the above results are somewnat encouraging, the requirement for 
refinements to the analysis is evident. These refinements should be incorpo- 
rated before the analysis is applied to determine airloadlngs for the impul- 
sive noise prollem. Recommendations regarding these refinements are presented 
below; 


1. Due to ccxnputer time and storsige limitations, application of the 
analysis using sufficiently small azimuth increments (e.g., 3 deg) to deter- 
mine the high harmonics of blade loading required for impulsive noise predic- 
tion is currently prohibitive. (The use of 15 deg increments for the ambient 
wind condition required approximately 3 hours of UNIVAC 1110 computer time.) 

A significant reduction in computer time and storage requirements could be 
achieved through selective vortex element interactions in the tip vortex 
geometry calculations such as described in Ref. 5. Also, the feasibility of 
acquiring higher harmonics of airloading without a general azimuth increment 
decrease should be explored. This might be achieved through locc.Uy 
decreasing the azimuth increment when a tip vortex is in close proximity to 

a blade, or through a decrease in the azimuth increment for the loading solu- 
tion relative to that used for the wake representation. 

2. In order to further reduce computer time and program complexity, a 
sensitivity study sb uld be performed to determine the cost effectiveness of 
several items, both input and internal to the computer program. For example, 
the sensitivity of the results to the input number of wake revolutions 
should be investigated to determine the proper compromise between accuracy 
and computer cost. Also, although the first flatwise bending mode deflec- 
tions contribute to the blade-vortex spacing, it is recognized that for many 
rotors the added complexity of including this aeroelastic blade response may 
be unwarranted. 
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3. The blade/wake representation should be refined in that it appears 
that the lifting line representation used may be insufficient for providing 
essential flow details at the blade tip. The sensitivity of the inflow at 
the tip tc the vortex representation in the inunediate vicinity of the blade 
should be investigated further. The influence of including a representation 
for shed vorticity and the sensitivity to time step size and stall limit 
should be included in the investigation. Regarding the requirement for a 
lifting surface representation of the blade, a better understanding of the 
I'undamental nature of impulsive noise would be helpful. That is, if it were 
known that impulsive noise cam be determined from chordwlse averaged air- 
loads, then lifting line assumptions may be sufficient. However, if it is 
dependent on local conditions on the blade, such as the occurrence of local 
pressure fluctuations or local sonic conditions, con^resslble lifting sur- 
face theory is required. 

4. Additional cases should be run to continue evaluation of the com- 
puter analysis. In particular, the analysis should be correlated with a 
test configviration for which a consistent set of wake geometrj' and airload 
data is acquired. Also, the analysis should be applied to transient gust 
conditions, to determine if transient blade-vortex Intersections arid asso- 
ciated impulsive type airloads are predicted which could be conducive to 
blade slap. 
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b) STRAIGHT LINE SEGMENT 


Geometry for Velocity Induced by a Vortex Segment. 
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Figure U. Scaling Factor to Approximate the Induced Velocity Contribution 
of an Adjacent Curved Vortex Segment. 
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a) END VIEW 



VORTEX CORE BOUNOARY 



b| SIDE VIEW 


Figure 5. 


Vortex Segment Geometry for Computing Velocities Within the 
Vortex Core. 
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HINGE OFF-SET 


Figure 6. Coordinate System. 
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i’igure 7. Gust Coordinates. 
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BLADf-WAKf REPflESENTATION 


REMARKS 
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Figure 9 - 


Blade-Wbke Vortex Sequent Representations . 
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Figure IG. Convergence of the Tip Vortex Geometry for a Hover Without Mind 
Condition and Comparison With Model Rotor Test Data. 
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BLADE THRUST LOADING, N/cm 
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(b) 0.939R and O .985 Blade Stations 


Figure 11. Concluded. 
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3 , (a) 0.75R and 0.85R Blade Stations. 
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z Figure 1?.. Final Calc\ilations of Time Histories of Blade Airloading for the Hovering Rotor With a 

i* Steady Aabient Wind Imposed at the First Time Step. 
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BLADE AZIMUTH ANGLE. DEG 

Fig^^^e 13. Time Histories of Blade Inplane Loading for the Hovering Rotor With a Steady Ambient Wind 
Imposed at the First Time Step. 
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Figure lU. Comparison of Predicted and Measured Blade-Tip Vortex Separation 
Distances for the Hovering Rotor With a Steady Ambient Hind. 
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